Correlation and isospin dynamics of participant-spectator matter in neutron-rich 

colliding nuclei at 50 MeV/nucleon 
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The sensitivities of isospin asymmetry and collision geometry dependences of participant (over- 
lapping region)-spectator (quasi-projectile and quasi-target region) matter towards the symmetry 
energy using the Isospin Quantum Molecular Dynamical model are explored. Particularly, the differ- 
ence of number of nucleons in overlapping zone to quasi projectile-target matter is found to be quite 
sensitive towards the symmetry energy at semi-peripheral geometries compared to the individual 
yield. It gives us a clue of this quantity to act as a measure for isospin migration. Further, the yield 
of neutrons (charge of the second largest fragment) are provided as a tool for overlapping region 
(quasi projectile-target) matter to check the sensitivity of above mentioned observable towards the 
symmetry energy experimentally. 

PACS numbers: 21.65.Ef, 25.70.Pq, 25.70.-z 



I. INTRODUCTION 

The primary goal of isospin physics is to obtain a better 
understanding of properties of asymmetric nuclear mat- 
ter or in other words isospin dependence of Nuclear Equa- 
tion of State (NEOS). Over three decades, the NEOS 
of symmetric nuclear matter is well understood by the 
study of giant dipole resonances, collective flow as well 
as multifragmcntation etc. [lH5j. The compressibility 

no = 9po d ^p A ^ Ip = po) which describes so called stiff- 
ness of the symmetric nuclear matter has been deter- 
mined to 235 ± 14 MeV. The NEOS of isospin asymmet- 
ric nuclear matter is recently underway, particularly, the 
density dependence of symmetry energy, which is crit- 
ically important also in many astrophysical processes, 
such as, structural and dynamical evolution of neutron 
star, the critical density for the direct cooling process 
etc. |3]. 

Considerable progress has been made in determining 
the sub- and supra-saturation density behavior of the 
symmetry energy 0-E3. The later part is still an unan- 
swered question in spite of recent findings in term of 
neutron-proton elliptic flow ratio and difference [HI, [13] . 
However, The former one is understood to some ex- 
tent [7Hll|. although, more efforts are needed for precise 
measurements. According to Symmetry Energy project 
(SEP) started from Michigan State University (MSU) at 
National Superconducting Cyclotron laboratory (NSCL) 
in 2011 (includes European as well as Asian countries), 
the sub-saturation density dependence of symmetry en- 
ergy still needs more study to verify the previously de- 
fined observables such as double ratio, isospin diffusion 
and then extend the study to define new observables to 
parametrize it (symmetry energy) more accurately, which 
has still long range in the literature from soft to linear 
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one @, H, Gil • Here we will try to search for new observ- 
able whose importance in intermediate energy heavy ion 
collisions is discussed below. 

Recently, the role of isospin degree of freedom has 
been investigated using collective flow and its balance en- 
ergy (at which flow disappears) [l6|. The collective flow 
is proven as an indicator for the symmetry energy At 
balance energy, the attractive interactions due to mean 
field are balanced by the repulsive interactions due to the 
nucleon-nucleon (NN) collisions. This counterbalancing 
is reflected in quantities like participant-spectator mat- 
ter [l7|. Recently, the participant-spectator matter at 
balance energy is found to be quite insensitive towards 
the mass and N/Z of the colliding system and hence 
can act as a barometer for the study of vanishing flow 
[l7l [l8j . The inter-dependence of collective flow with 
symmetry energy and collective flow with participant- 
spectator matter gives us a clue to check the sensitivity of 
participant-spectator matter towards the symmetry en- 
ergy in asymmetric colliding nuclei. 

Secondly, the elliptic flow is also shaped by the inter- 
play of collisions and mean field [l9| . In addition, the el- 
liptic flow pattern of the participant matter is affected by 
the presence of cold spectator matter (20[ ■ Specially, the 
spectator can inhibit the collective transverse expansion 
of the decompressing participant matter and effectively 
shadow particle emission directed towards the reaction 
plane. This study is indicating that elliptical flow is influ- 
enced strongly by participant-spectator matter distribu- 
tions. As we have discussed earlier, the elliptic flow ratio 
and difference from neutrons and protons was used as an 
indicator for the symmetry energy [ill [l4| . It is the an- 
other insight/clue to check the sensitivity of participant- 
spectator matter towards the symmetry energy. 

Apart from collective transverse and elliptic flow, the 
participant-spectator matter also plays an important role 
in understanding multifragmentation as well as nuclear 
stopping. In recent years, the correlation between nuclear 
stopping and light charged particles (LCPs) is investi- 
gated by using Quantum Molecular Dynamics (QMD) 
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and Isospin-QMD (IQMD) [21|. The relation is estab- 
lished between nuclear stopping, directed flow and elliptic 
flow [2^| . Further, the participant matter is declared as 
an indicator to nuclear stopping [l7| . All the correlation 
studies indicate some indirect correlation of participant- 
spectator matter with different kind of fragments. So, it 
also becomes our prime duty to correlate the participant- 
spectator matter with different kind of fragments and 
then provide specific type of fragments as an experimen- 
tal measure for participant-spectator matter. 

In the literature, mostly participant-spectator termi- 
nology is used, which is most suitable at higher inci- 
dent energies, where the dynamics is almost fully ac- 
counted by NN collisions. In the energy region of present 
manuscript, the reaction proceeds through the interplay 
of both NN collisions and mean field. We prefer to use 
quasi-projectile (QP) or quasi-target (QT) for spectator 
matter and overlap region (OR) for the participant mat- 
ter. 

In the present work, IQMD model is used, which is 
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discussed in detail in our recent publications 23 1, 
originally developed by Hartnack and collaborators 24 1 . 
The model is modified by the authors for the density 
dependence of symmetry energy, having form: 



Esym(p) 



C s 



2/3 



c s 



with the parameters of C Sj fe = 25 MeV and C SiP = 35.2 
MeV. When we set 7^ = 0.5 and 1.5, respectively, it 
corresponds to the soft and stiff symmetry energy 



II. RESULTS AND DISCUSSION 

The present OR and QP+QT matter demonstration is 
based on the fireball model [2f| as reported in Ref . fl7l] . 
All nucleons having experienced at least one collision are 
supposed to originate from OR matter (labeled as Nor)- 
The remaining matter is called QP,QT matter (labeled 
as Nqp+qt). The Nxot are the total number of nucleons 
in the reaction system. This concept gives similar results 
as has been demonstrated in the fireball model of Gosset 



et al. [25[ and further verified by QMD in last couple 
of years to till date 0, EH- There is another way to 
define the OR, QP and QT matter on the basis of differ- 
ent rapidity cuts @,[l7|- However, both these definitions 
have been reported to give the same results d, [l7| ■ In 
the present study, we shall use the first definition to con- 
struct the OR and QP+QT matter. During the expan- 
sion stage, this definition will lead to the production of 
two matters with different densities. The OR must have 
relatively low density compared to QP,QT matter. This 
density gradient must increase, when one goes towards 
the peripheral collisions. Due to the density gradient, 
the transfer of particles from high to low density region 
is found to relate with the phenomenon known as isospin 
migration [26{ . 
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FIG. 1: (Color online) Impact parameter (left panels) and 
isospin asymmetry (right panels) dependences of nucleons in 
OR (a,b) and QP+QT (c,d) matter including the contribu- 
tions from neutrons and protons. The solid (open) symbols 
represent soft (stiff) symmetry energies. The OR, QP and QT 
are Overlapping Region, Quasi-Projectile and Quasi- Target, 
respectively. 



Mathematically, isospin migration can be understood 



as: DP -DP cx id 



dE s 



withC, the mass coefficients, 

p/n 



which are directly given by the variation of n,p chemical 
potentials with respect to density and asymmetry. For 
more detail see Ref. (2(|. From the above equation, it 
is clear that isospin migration depends on the slope of 
the symmetry energy, or the symmetry pressure. When 
one moves towards the semi-peripheral collisions, QP or 
QT of about normal density are in contact with the OR 
where the density is quite low than saturation density. 
In this region of density, a stiff symmetry energy has a 
smaller value but a larger slope in comparison with a soft 
symmetry energy. In brief, this definition can help us to 
find the observable for isospin migration. 



A. Theoretical probe for isospin migration 

In the present study, thousand of events are simulated 
using static soft equation of state and energy dependent 
NN cross-section for the isotopes of Sn, namely 112 Sn + 
112 Sn, 124 Sn + 124 Sn and 132 Sn + 132 Sn at incident en- 
ergy 50 MeV/nucleon along the whole collision geometry. 

Fig. [1] displays the impact parameter (b) and isospin 
asymmetry of the system ((N/Z) system ) dependences of 
OR and QP+QT including the neutrons and protons con- 
tribution. With the increase of impact parameter, the 
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number of nucleons in OR (QP+QT) matter decreases 
(increases). The similar behavior is followed by the im- 
pact parameter dependence of neutrons and protons. In 
contrast, there exists more sensitivity for the number dif- 
ference between the neutrons and protons near the cen- 
tral (peripheral) geometries from OR (QP+QT) matter, 
which is due to the dominance of more NN collisions (very 
less collisions) near the respective collision geometries. 

In semi-peripheral collisions, the isospin asymmetry 
dependence of OR and QP+QT matter along with neu- 
trons and protons contribution have dramatic behavior. 
With the increase in isospin asymmetry, the number of 
nucleons in OR (QP+QT) matter increases (decreases). 
Furthermore, the neutrons (protons) contribution from 
OR as well as QP+QT matter increases (decreases). In- 
terestingly, the sensitivity is stronger for QP+QT matter 
at higher isospin asymmetry. This is due to the fact that 
(1) With isospin asymmetry, there emerges a sharp in- 
crease in neutron content inside QP+QT about 4-5 times 
compared to the decrease of QP+QT matter (see the red 
circles vs black circles in right bottom panel); (2) As the 
proton content is constant for all the isotopes of Sn and 
while the QP+QT matter is decreasing, it will eventually 
lead to the decrease in protons content in the QP+QT 
matter. 

Lastly, the effect of symmetry energy is relatively weak 
on the impact parameter and isospin asymmetry depen- 
dences on the nucleons originated from OR and QP+QT 
matter. This difference is true as the error bars are less 
than the size of the symbols. The soft (stiff) symmetry 
energy contributes more for OR (QP+QT) matter. This 
is due to the gradient in densities of two matters, which 
can be explained by slope of symmetry energy rather 
than magnitude. The slope (magnitude) of symmetry 
energy changes its behavior below (at) the saturation 
density for the soft and stiff symmetry energy. Due to 
the large slope (less magnitude) for the stiff symmetry 
energy even below the saturation density, the QP+QT is 
more neutron-rich with the stiff symmetry energy, which 
is true with the soft symmetry energy for OR matter as 
the density at freeze-out time in comparison to QP+QT 
is very low. 

From the above discussions, it is clear that number of 
nucleons contribution from OR and QP+QT matter is 
a good candidate to explore the isospin physics, but not 
a potential candidate for symmetry energy (due to weak 
dependence). However, the density gradient of symmetry 
energy gives us some preliminary clues of isospin migra- 
tion between the nucleons of OR and QP+QT matter. 

In order to reveal the effect of symmetry energy in 
term of isospin migration, in Fig. [21 we display the 
difference of nuclcon number between OR and QP+QT 
( Af °' R ^ f F+QT ) as a function of kinetic energy of nucleons 
at different impact parameter and for two neutron-rich 
systems. In central collisions, no effect of symmetry en- 
ergy is observed. However, with the increase of impact 
parameter, sensitivity of symmetry energy towards the 
observable seems increasing. This indicates that the ef- 
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FIG. 2: (Color online) The kinetic energy spectra of the 
Nor JVqp+qt different colliding geometries. The solid 
(dashed) lines represent the contribution with soft and stiff 
symmetry energy. The a, b, c, d and e, f, g, h are for 
112 Sn+ n ' 2 Sn and 124 Sn+ 124 Sn systems, respectively. 

feet of symmetry energy on the observable under study 
depends weakly on the isospin asymmetry, but strongly 
on the density gradient. As discussed earlier, the density 
gradient increases strongly (weakly) with the impact pa- 
rameter (isospin asymmetry). The effect of symmetry en- 
ergy on the observable ( NoR ~ 1 ^ r Q ^ +QT ) mainly originates 
from the density difference and not due to the isospin 
asymmetry. We further know, the density difference is a 
direct measure of the isospin migration. The findings are 
also supporting the results of the Ref. (26[. 

From the study, we can conclude that the difference 
of number of nucleons from OR to QP+QT can act as a 
probe to the slope of symmetry energy versus the density 
in term of isospin migration at sub-saturation densities. 



B. Correlation with fragments and observables for 
experiments 

The second aspect of this paper is to correlate the num- 
ber of nucleons in OR and QP+QT matter with fragmen- 
tation process and then provide some observables to test 
the above observable as a sensitive probe for symmetry 
energy. 

To this end, we displayed the impact parameter and 
isospin asymmetry dependences of different kind of frag- 
ments in Fig. [31 which is just like in Fig. [T] The 
multiplicity of different kind of light charged particles 
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FIG. 3: (Color online) Same as Fig. [1] but for the multiplic- 
ity of different kind of fragments from lighter to the heavier 
ones. In upper panels (a,b), the solid squares represent for 
free nucleons, solid circles (open squares) for free neutrons 
(protons), circles (triangles) for deuteron (triton), open stars 
(solid stars) for 3 He ( 4 He). In lower panels (c, d), solid and 
open diamonds for the heaviest fragment and second heaviest 
fragment, respectively. 



goes on decreasing with the increase of impact param- 
eter as well as size of the fragments. This type of the 
behavior has been already observed many times in the 
literatures 
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On the other hand, Z max _i/A pro j ectae 
(charge of the second largest fragment normalized by the 
mass number of projectile) and Z max / A pro j ect ii e (charge 
of the largest fragment normalized by the mass number 
of projectile) increase with impact parameter. These two 
L. H.S. panels are behaving exactly the similar just like 
the L. H.S. panels in Fig. [TJ indicating some clues of 
the correlation for the measure of number of nucleons in 
OR(QP+QT) using the light charged particles (charge of 
the heavier fragments). 

To further strengthen the correlation, the R.H.S. 
of Fig. [3] give us some interesting features. Just 
like the isospin asymmetry dependence of nucleons in 
OR (QP+QT) matter, the multiplicity of free nucleons 
(charge of the heavier fragments) is increasing (decreas- 
ing). Moreover, with the isospin asymmetry, the increas- 
ing (decreasing) trend of neutrons (protons) contribution 
in OR matter (FigJT]) is similar to the behavior of the mul- 
tiplicity of neutrons (protons) (Fig|3]). In addition, the 
isospin asymmetry dependence of charges of heavier frag- 
ments is also similar to the contribution of protons from 
the QP+QT (Fig. [T]). These are the strong evidences 
for correlation. The important point to check here is 
the sensitivity of yield of different kind of fragments and 
charge of the heavier fragments with nucleons of OR and 
QP+QT matter, which will be further checked in Fig. @] 
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FIG. 4: (Color online) Correlation between the nucleons of 
OR and QT+QT matter with different kind of fragments. 
The lines are fitted with the power-law form Y = CX T . Sym- 
bols are the same as Fig. 3. 



TABLE I: r value extracted from Fig(?l showing the positive 
(negative) slope of light charged particles with participant 
(spectator) matter. 



Particles 


n 


P 


' 2 R 


3 H 


* He 


4 He 


t for OR 
t for QP+QT 


0.605 
-0.317 


0.546 
-0.288 


0.815 
-0.423 


0.831 
-0.418 


1.022 
-0.524 


1.180 
-0.599 



by using the power-law fitting method. 

From the light charged particles, interesting isospin ef- 
fects arc observed with isospin asymmetry of the sys- 
tem. With the increase in the one neutron in the LCPs, 
i.e. from 1 H— >- 2 H— » 3 H and 3 He— ^ 4 He, the multiplicity 
changes its trend from decreasing 1 H, 3 He to increas- 
ing for 3 H, 4 He, respectively. However, with the in- 
crease of one proton from 3 H to 3 He, the multiplicity 
reveals a sharp change from increasing to decreasing. 
The neutrons-protons distributions within the LCPs is 
also satisfying the neutrons-protons distributions from 
the OR zone. This indicates the strong correlation be- 
tween LCPs and OR matter, but the point of interest is 
which type of LCPs are perfect indicator for OR matter 
which will be discussed in FigPJJ 

In Fig. 21 the sensitivity of number of nucleons in OR 
and QP+QT matter with the multiplicity of LCPs as 
well as charge of the heavier fragments is checked by fit- 
ting the power law of the form Y = CX T . The positive 
correlation is observed between OR (QP+QT) nucleons 
and multiplicity of LCPs (Z max and Z max -i). Interest- 
ingly, from the power-law slope, it is found that although, 
the multiplicity decreases with the size of light fragments, 
the slope parameter or sensitivity increases with the frag- 
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ment size (shown in Tabic [TJ) . This type of prediction is 
also true for the second largest fragment and the largest 
fragment for the QP+QT nucleons. The slope parame- 
ter for neutrons (0.6) with nucleons of OR matter (also 
shown in Table [XJ) and for Z max -i/Ap ro j ect u e with nu- 
cleons of QP+QT matter (0.58) (FiglJ) is almost the 
same, which reveals the similar sensitivity of the respec- 
tive matter towards the respective fragments. From here, 
one can say that if one uses the neutron as a measure for 
the nucleons of OR matter and Z max _i/ Ap ro j ecti i e as a 
nucleons of QP+QT matter, then the difference of two 
parameters at semi-peripheral geometry, just like in Fig. 
[2l can act as a probe for the isospin migration or slope 
of symmetry energy. From the Table HI we also observed 
that there is a linear correlation between the nucleons of 
OR matter and 3 He fragments (i.e., having r very close 
to one), which indicates that the 3 He can be taken as a 
direct measure for the nucleons of OR matter. 



III. CONCLUSION 

In summary, we tried to find an observable to measure 
density dependence of symmetry energy or its slope by 



studying the OR and QP+QT matter in intermediate en- 
ergy heavy-ion collisions. The difference in number of nu- 
cleons of OR to QP+QT matter is particularly sensitive 
towards the slope of symmetry energy at semi-peripheral 
geometries. This gives us a clue that this observable can 
act as probe for the isospin migration. Principally, the 
yield of neutrons and charge of the second largest frag- 
ment {Z max -i/Ap ro j ect u e ) could be provided as experi- 
mental tools to check the sensitivity of nucleons in OR 
and QP+QT matter towards the symmetry energy in 
term of isospin migration. Interestingly, the nucleons of 
OR matter has a linear correlation with the yield of 3 He 
(t close to 1). 
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